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I. INTRODUCTION 

The theoretical understanding of neutron star (NS) struc- 
ture requires the knowledge of the equation of state (EOS) 
of highly compressed cold baryonic matter, up to densities of 
about ten times normal nuclear density, po ~ 0.17 fm~^ iH. 
In such an extreme environment, the appearance of "exotic" 
components of matter, such as hyperons, meson condensates, 
and quark matter (QM), is expected |2]. 

It is in fact well known that hyperons appear at around 
2-3 po in beta-stable nuclear matter and lead to a strong soft- 
ening of the EOS with a consequent substantial reduction of 
the maximum NS mass Ist |4j]- The theoretical maximum 
mass of hyperonic NS can even result below current obser- 
vational values of about 1.5 solar masses jst], as in the case 
of the microscopic Brueckner-Hartree-Fock (BHF) approach 
for the hyperonic EOS [6, JJ, which we employ in this work 
for the hadronic phase. This would imply the presence of 
nonhadronic "quark" matter in the interior of massive NS 
is, 9, 10, 1 1], heavier than the maximum mass hyperon stars. 

However, the appearance of quark matter poses the prob- 
lem of an accurate theoretical description of the quark phase, 
which is so far an open question, and furthermore of the de- 
tails of the phase transition between hadronic and quark mat- 
ter. The purpose of this article is the study of the latter prob- 
lem, combining the BHF EOS of hyperonic matter with a 
generalized phenomenological MIT bag model for the quark 
phase. In our case the EOS of hyperonic matter is so soft that 
the deconfinement transition occurs at rather low densities, 
where hyperon contamination is not so large. Assuming the 
quark deconfinement transition to be of first order, it causes 
a thermodynamic al instability and the mixed phase (MP) ap- 
pears around the critical density. In the usual Maxwell con- 
struction (MC), the MP is composed of two charge-neutral 
hadron and quark phases and uniform density distributions are 
assumed in each phase. 

The properties of the MP are very different, if there are 
more than one independent chemical potentials, as in our case. 
Glendenning has pointed out that in this case the MC is not the 
coiTect procedure to obtain a thermodynamically well-defined 
EOS with the MP, but that one must properly satisfy more 
fundamental requirements by means of the Gibbs conditions 



III2II . Since then many works have appeared regarding nu- 
clear pasta structures in low-density nuclear matter L13„ .I4l1 , 
kaon condensation at several times po ITs^ 14 , 3, IBl' ^"^d the 
hadron-quark deconfinement transition 
When the Gibbs conditions are applied to the quark deconfine- 
ment transition, the MP is composed of individually charged 
hadron and quark phases, and baryon number density as well 
as charge density are nonuniform in each phase, arranged in 
different geometrical structures. 

However, electromagnetic and surface contributions to the 
energy of the MP are only approximately treated in the usual 
bulk calculations [VT}, but could have an important effect 
lITil [161 [Tsi [Toll . The quantitative analysis of these correc- 
tions (sometimes called finite-size effects) for the quark de- 
confinement transition in hyperonic matter is one purpose of 
this article. We elaborately figure out the roles of the finite- 
size effects in the mixed phase. We shall see that these effects 
change remarkably the properties of the mixed phase; e.g., the 
geometrical structures are destabilized by the charge screen- 
ing effect for the Coulomb interaction in the extreme case, and 
the EOS given by the MC is effectively recovered. Regarding 
hyperon mixing we shall see that the appearance of the mixed 
phase completely suppresses the appearance of hyperons. 

Some major results of our work were already presented in 
another paper |20], while here we provide more details and 
furthermore study the influence of the mixed phase on the 
global NS observables like mass-radius relation and in par- 
ticular the maximum mass. We shall see that the global prop- 
erties of compact stars are little changed, but the structure and 
property of the internal core are very different, compared to 
the MC. Such difference may affect the elementary processes 
like neutrino transport or neutrino emission in the core. 

In the following, we give in Sec. a concise summary of 
the BHF approach for hyperonic matter that is used, whereas 
in Sec. |III]the modified MIT bag model for the quark phase is 
introduced. Sec. |IV] discusses the details of the properties of 
the mixed phase and figures out the peculiar role of the finite- 
size effects. In Sec. [V] the EOS including the quark decon- 
finement transition is applied to the structure of hybrid stars. 
Summary and concluding remarks are given in Sec. |Vll 
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n. BHF APPROACH OF HYPERONIC MATTER 

Our theoretical framework for the hadronic phase of mat- 
ter is the nonrelativistic Brueckner-Hartree-Fock approach 
ll2lll based on microscopic nucleon-nucleon (NN), nucleon- 
hyperon (NY), and hyperon-hyperon (YY) potentials that are 
fitted to scattering phase shifts, where possible. Nucleonic 
three-body forces (TBF) are included in order to (slightly) 
shift the saturation point of purely nucleonic matter to the em- 
pirical value. 

It has been demonstrated that the theoretical basis of the 
BHF method, the hole-line expansion, is well founded: the nu- 
clear EOS can be calculated with good accuracy in the BHF 
two hole-line approximation with the continuous choice for 
the single-particle potential, since the results in this scheme 
are quite close to the full convergent calculations which in- 
clude also the three hole-line contributions |22]. Due to these 
facts, combined with the absence of adjustable parameters, the 
BHF model is a reliable and well-controlled theoretical ap- 
proach for the study of dense baryonic matter 

In the following we give a short review of the BHF ap- 
proach including hyperons. Detailed accounts can be found 
in Refs. S|23l|23l. The basic input quantities in the Bethe- 
Goldstone equation are the NN, NY, and YY potentials. In 
this work we use the Argonne Vis NN potential [2^] supple- 
mented by the Urbana UIX nucleonic TBF of Ref . i26tl and the 
Nijmegen soft-core NSC89 NY potentials fl^l that are well 
adapted to the existing experimental NY scattering data and 
also compatible with A hypernuclear levels ll28l[29tl . Unfortu- 
nately, the NSC89 potentials contain no YY components, be- 
cause up to date no YY scattering data exist. Nevertheless the 
importance of YY potentials should be minor as long as the 
hyperonic partial densities remain limited. Recently also cal- 
culations using the NSC97 NY and YY potentials Q were 
completed, which yield very similar maximum NS masses in 
spite of quite different internal compositions |7]. 

With these potentials, the various G matrices are evaluated 
by solving numerically the Bethe-Goldstone equation, which 
can be written in operator form as 

Qc 



Gab[W] = Vab 



(1) 

where the indices a, 6, c indicate pairs of baryons and the Pauli 
operator Q and energy E determine the propagation of in- 
termediate baryon pairs. The pair energy in a given channel 
c = = A,S- is 



(2) 



with Ti{k) = riii + fc^/2mi, where the various single-particle 
potentials are given by 



(3) 



J— n,p.A,S~ 

and are determined self-consistently from the G matrices, 
Ul'\k) = ^ Re(fcfc'|G(,,)(y) k')] \kk') . (4) 



The coupled equations ([T]i to (|4| define the BHF scheme with 
the continuous choice of the single-particle energies. In con- 
trast to the standard purely nucleonic calculation, the addi- 
tional coupled-channel structure renders the calculations quite 
time-consuming. 

Once the different single-particle potentials are known, the 
total nonrelativistic hadronic energy density, e^f, can be eval- 
uated: 



E 



1, 



(5) 



and €h is thus represented as a function of particle number 
densities piii — n,p, A, S^) for a given baryon number den- 
sity pb- Knowing the hadronic energy density, and adding 
the contributions of the noninteracting leptons (l = e, /i), 
e = ch + the various chemical potentials pi = de/dpi (of 
the species i — n,p,A,Y,~,e,p) can be computed straight- 
forwardly and the equations for beta-equilibrium. 



Pi B^Pn QiPe 5 



(6) 



{Bi and Qi denoting baryon number and electric charge of 
species i), baryon number density and charge neutrality. 



'^BiPi = Pb 

i 

^QtPt = 0, 



(7) 
(8) 



allow one to determine the equilibrium composition Pi{pB) 
at a given baryon number density pB and finally the hadronic 
EOS, 

p ^ N 2 d e{p,{pB)) de 

PH{pB)^PB-i = PB-] e- (9) 

dpB Pb apB 



III. QUARK MATTER EOS 

Unfortunately, the current theoretical description of quark 
matter is burdened with large uncertainties, seriously limiting 
the predictive power of any theoretical approach at high bary- 
onic density. For the time being we can therefore only resort 
to phenomenological models for the quark matter EOS and try 
to constrain them as well as possible by the few experimental 
information available on high-density baryonic matter. 

One important condition is due to the fact that certainly in 
symmetric nuclear matter no phase transition is observed be- 
low w 3/3o- We will in the following use an extended MIT bag 
model |31] (requiring a density-dependent bag "constant") 
that is compatible with this condition. 



A. The MIT bag model 

We first review briefly the description of the bulk proper- 
ties of uniform quark matter, deconfined from the beta-stable 
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hadronic matter mentioned in the previous section, by using 
the MIT bag model [31i] . The energy density of the f = u,d,s 
quark gas can be expressed as a sum of the kinetic term and the 
leading-order one-gluon-exchange term ll32ll33ll for the inter- 
action energy proportional to the QCD fine structure constant 



ical potentials are modified as 



eg = 



5 + 

/ 

Xf (2a; J + l) yjl + x'j — arsinhx/ 



(10) 



3my 
8^ 



'1 + xj 



rsinh 



(11) 



where m/ is the / current quark mass, Xf = kl^P/mf, the 

number density of / quarks is p/ = k^p^ / tt^, and B is the en- 
ergy density difference between the perturbative vacuum and 
the true vacuum, i.e., the bag "constant." We use massless u 
and d quarks and — 150 MeV. 

It has been found lH, IM] that within the MIT bag model 
(without color superconductivity) with a density-independent 
bag constant B, the maximum mass of a NS cannot exceed 
a value of about 1.6 solar masses. Indeed, the maximum 
mass increases as the value of B decreases, but too small 
values of B are incompatible with a hadron-quark transition 
density pB > 2-3 po in nearly symmetric nuclear matter, as 
demanded by heavy-ion collision phenomenology. Values of 
B > 150 MeV/fm'^ can also be excluded within our model, 
since we do not obtain any more a phase transition in beta- 
stable matter in combination with our hadronic EOS |8]. 

In order to overcome these restrictions of the model, one 
can introduce empirically a density-dependent bag parame- 
ter Bcffips), which has not any more the simple interpre- 
tation as the energy difference between the perturbative and 
the true vacua; instead some density dependence originat- 
ing from the non-perturbative interaction energy AeintiPs) 
is considered to be included in the effective bag constant, i.e., 
-Scff(ps) — B + Aeint(pB) 1 10]. This allows one to lower 
the value of B at large density, providing a stiffer QM EOS 
and increasing the value of the maximum NS mass, while at 
the same time still fulfilling the condition of no phase transi- 
tion below PB ~ 3/30 in symmetric matter. In the following 
we present results based on the MIT model using both a con- 
stant value of the bag parameter, B = 100 McV/fm'^, and a 
Gaussian parametrization for the density dependence. 



Bcff(pB) = S+(Bo-B)cxp 



^ Po ^ 



(12) 



with B = 50 MeV/fm3, Bq = 400 McV/fm^, and /? = 0.17. 
For a more extensive discussion of this topic, the reader is 
refen-ed to Refs. jslUfToll. 

The introduction of a density-dependent bag parameter has 
to be taken into account properly for the computation of vari- 
ous thermodynamical quantities; in particular the quark chem- 



1 dBcsipB) 



dpi 



(13) 



Nevertheless, due to a cancellation of the second term in ( fT3] ), 
occurring in relations (fT4l l for the beta-equilibrium, the com- 
position of QM at a given total baryon number density remains 
unaffected by this term (and is in fact independent of B). At 
this stage of investigation, we disregard possible dependencies 
of the bag parameter on the individual quark densities. 

In the beta-stable pure quark phase, the individual quark 
chemical potentials are fixed by Eq. (|6| with Bf = 1/3, which 
implies 



(14) 



The charge neutrality condition and the total baryon number 
conservation read 

2 1 1 

= -PU - - - Pi , (15) 

1 . 

PB = g [Pu + Pd + Ps) ■ (16) 

These equations determine the composition Pf{pB) and the 
pressure of the QM phase. 



Pq^Pb) = Pb 



2 dieq/pB) 



(17) 



The modified bag model is clearly an oversimplified model 
of QM, but in this article we focus mainly on the differences 
between the different hadron-quark phase transition construc- 
tions in NS matter introduced in the following. 



IV. HADRON-QUARK PHASE TRANSITION 

A. Gibbs conditions and the mixed phase 

Figure [T] compares the hadronic BHF EOS and the quark 
matter EOS with different values of the parameters B and 
Q!s for beta-stable and charge-neutral matter. One can see 
that the quark EOS approaches that of a relativistic free gas 

1/3 

(E/A ^ p^ ) with increasing density, while the hyperonic 
EOS is always soft. Consequently the quark deconfinement 
transition cannot occur at too high densities. If we demand the 
quark and the hyperonic EOS to cross, ag should be small and 
B slightly large, which gives a relatively low critical density. 
Thus the appearance of hyperons is effectively suppressed due 
to a quark deconfinement transition. In this article we choose 
as = and compare results with B = 100 MeV/im^ and 
Bcs{pb)- 

Since we shall see that two independent chemical poten- 
tials, charge chemical potential and baryon number chemical 
potential, are needed to specify the thermodynamical proper- 
ties of the MP, we must properly take into account the Gibbs 
conditions |17], which require the pressure balance and the 
equality of the chemical potentials between the two phases 
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FIG. 1: (Color online) EOS of beta-stable hadronic matter (black 
curves) and quark matter (colored curves) with constant B — 
60 MeV/fm'^ (lower curves) and B = 100 MeV/fm"^ (upper curves) 
for several values of Us. The gray curve shows the result for the 
Bc«{pb) model with = 0. Hyperons appear at the dotted point 
in hadronic matter. 



besides the thermal equilibrium. We employ a Wigner-Seitz 
approximation in which the whole space is divided into equiv- 
alent cells with a given geometrical symmetry, specified by the 
dimensionality d = 3 (sphere), d — 2 (rod), or d ~ 1 (slab). 
In each cell with volume Vw the quark and hadron phases are 
separated: a lump portion made of the quark phase with vol- 
ume Vq is embedded in the hadronic phase with volume Vh 
or vice versa. A sharp boundary is assumed between the two 
phases and the surface energy is taken into account in terms 
of a surface-tension parameter <j. The surface tension of the 
hadron-quark interface is poorly known, but some theoretical 
estimates based on the MIT bag model for strangelets i33[l and 
lattice gauge simulations at finite temperature |35t] suggest a 
range of cr « 10 — 100 MeV/fm^. We show results using 
cr = 40 McV/ fm^ in this article, and discuss the effects of its 
variation. 

We use the Thomas-Fermi approximation for the density 
profiles of hadrons and quarks. The energy for each cell is 
then given as 

E= ( SreH[p^{r)] + ( d^r tQ[pf{r)] + E, + Ec + aS 

JVh JVq 

(18) 

with i — n,p, A, E^, f = u,d, s, and S being the hadron- 
quark interface area. E^ is the electron gas kinetic energy. 



JVw 



(19) 



approximated by ee[Pe] — {^tt'^ Pe)^^^/ 4:n^ . For simplicity 
muon is not included in this article. Ec is the Coulomb inter- 
action energy, 

I' d'^rd^'r' ^-h(^)Pch(r') 



Er 



Vw 



\r — r' 



(20) 



where the charge density is given by epch(f") = 
T,^=n,p,A,^-,eQ^P^ir) m Vh and epch(r) 



J2f=uA.s,eQfPfi'^) ™ with Qa being the particle 
charge (Q = —e < for the electron). Accordingly, the 
Coulomb potential 0(r) is defined as 



</>(r) 



^3^/ e^Pchjr') 



(21) 



Vw 



where 0o is an arbitrary constant representing the gauge de- 
gree of freedom. We fix it bystipulating the condition 
(f>{Rw) — 0, as before lfl3l ITsl TToll. The Poisson equation 
then reads 



A0(r) = 47reVch(r) 



(22) 



The density profiles of the hadrons, Pi{r), i = n,p, A, S", 
and quarks, Pf{r), f = u,d, s, are determined by the equa- 
tions of motion, 

S{E/Vw) 



Spa{r) 



Pa , 



(23) 



where we introduced the chemical potentials p.a for the par- 
ticle species a — n,p,A,Y,^,u,d,s,e. Note that some ad- 
ditional terms are needed for the quark chemical potentials 
in the case of Bcs{pb) [see Eq. (fTSTll. We consider chemi- 
cal equilibrium at the hadron-quark interface as well as inside 
each phase, so that Eq. ^ implies 

p.u+ Pe= Pd = Ps , 

Pp + Pe = Pn — PA — Pu + 2pd , 

Ps- + Pp = "^Pn ■ (24) 

There are two independent chemical potentials, which are usu- 
ally chosen as the charge chemical potential pq = p^ and the 
baryon number chemical potential pB — Pn- For a given 
baryon number density 



PB 



1 



E 

i— n,p,A.S 



d^rp,{r) + 



E 

f—u,d,s 



d^' 

Vq 



.Pfjr) 



(25) 

Eqs. (|221424] |. together with the global charge neutrality con- 
dition, Jy^^ d^rpch{r) = 0, obviously fulfill the requirements 
by the Gibbs conditions. The optimum dimensionality d of 
the cell or the lump, the cell size Rw, and the lump size R 
[or equivalently the volume fraction fy = {R/ Rw)"^]^ are 
searched for to give the minimum energy. 

Note that the Poisson equation (l22l i is a highly non- 
linear equation with respect to the Coulomb potential through 
Eq. ( |23] |. If we linearize it, we obtain the Debye screening 
length 



dpi 



(26) 



(with i = n, p, A, , e and i = u,d,s,e in the hadron and 
quark phases, respectively), which gives a rough measure to 
estimate how effective is the charge screening. We shall see 
that Ad ^ ©(several fm) < R, Rw, which confirms the im- 
portance of the screening effect. We use the relaxation method 
to solve Eqs. ( |221 - |24] | consistently. The details of the numeri- 
cal procedure to calculate the EOS and determine the geomet- 
rical structure of the MP are similar to that in Refs. lfl3l[T9ll . 
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FIG. 2: (Color online) Upper panel: Density profiles and Coulomb 
potential (ji within a 3D (quark droplet) Wigner-Seitz cell of the MP 
at pb ~ 0.4 fm""^. The cell radius and the droplet radius are Rw = 
21.6 fm and R — 13.8 fm, respectively. Middle panel: Density 
profiles in the MC case, where only the volume fraction is calculated, 
while the absolute values of R and Rw are irrelevant. We use the 
same Rw as in the upper panel, while the volume fraction of the 
quark phase is different, fv = 0.375. See Fig. [3] for the volume 
fraction. The results are for the B = 100 MeV/fm^, as = 0, 
a — 40 MeV/fm^ case. Lower panel: Density profiles in the bulk 
Gibbs calculation. 



B. Charge screening effect and the Maxwell construction 

Figure |2] shows an example of the density profile in a 3D 
cell for Pb = 0.4 fm"^ and the B = 100 McV/fm^ case. 
One observes a non-uniform density distribution of each par- 
ticle species together with the finite Coulomb potential, while 
bulk calculations use a uniform density without the Coulomb 
interaction Jsl [13, [13] • This is due to the charge screening 
effect: the charged particle distribution is rearranged to give 




Pe[fm 



FIG. 3: (Color online) Wigner-Seitz cell radius Rw (upper colored 
segments), droplet radius R (lower colored segments), and volume 
fraction fv (black segments) as a function of baryon density. The 
results are for the B = 100 MeV/fm^, = 0, cr = 40 MeV/fm^ 



smaller Coulomb interaction energy. Different from the MC 
case (middle panel), where local charge neutrality is implic- 
itly assumed, the quark phase is negatively charged, so that 
d and s quarks are repelled to the phase boundary, while u 
quarks gather at the center. Thus local charge neutrality is ob- 
viously violated at any point inside the cell, even at the center 
of the droplet or at the boundary of the cell. The hadron phase 
is positively charged: protons are attracted by the negatively 
charged quark lump, while electrons are repelled. We shall 
see that such reaiTangement gives rise to a remarkable effect 
on the energy of the MP. 

The density dependence of the optimal structures and their 
characteristic dimensions R, Rw, and the volume fraction 
fv are shown in Fig. [3] One observes a transition from 
droplet to slab to tube to bubble with increasing density. With 
B = 100 MeV/fm^ and Bcs{pb), the MP occurs in the inter- 
val PB = 0.298 - 0.708 and pB = 0.236 - 0.670 fm^^ 
respectively, i.e., it appears at less than twice normal nuclear 
density and extends up to much larger density, relevant for 
NS physics. The transitions between the different geometric 
structures are by construction discontinuous and a more so- 
phisticated approach would be required for a more realistic 
description of this feature. 

In the lower panel of Fig.|2] shown is the case of bulk Gibbs 
calculation. The local charge neutrality is notreaUzed as in the 
full calculation. However, by ignoring the Coulomb interac- 
tion, the density distribution in each phase is uniform. 

Figure ID (upper panel) compares for B = 100 MeV/fm"^ 
the resulting energy per baryon of the hadron-quark MP with 
that of the pure hadron and quark phases over the relevant 
range of baryon density. The thick black curve indicates the 
case of the MC, while the colored line indicates the MP in its 
various geometric realizations by the full calculation. While 
the structure and the composition of the MP by the full calcu- 
lation are very different from those of the MC case, as shown 
before, the energy of the MP is only slightly lower than that 
of the MC, and the resulting EOS is similar to the MC one. 

The pressure is presented in the lower panel of Fig. ID In 
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FIG. 4: (Color online) EOS of the MP (thick curves) in compar- 
ison with pure hadron and quark phases (thin curves). The upper 
panel shows the energy per baryon EjA, the middle panel the en- 
ergy difference between mixed and hadron (pB < 0.44 fm~^) or 
quark (ps > 0.44 fm~'^) phases, and the lower panel the pressure. 
Different colored segments of the MP are chosen by minimizing the 
energy. The EOS within the MC (between pn = 0.34 fm^'' and 
pQ — 0.50 fm~"^) is also shown for comparison (thick black curves). 



the MC case, where local charge neutrality is implicitly as- 
sumed, the two equilibrium densities pn ~ 0.34 fm^"^ and 
Pq w 0.50 fm^-^ of the pure phases are connected by a 
straight line with equal pressure Ph{ph) = Pq{pq)- Com- 
pared with the MC, the pressure of the MP smoothly inter- 
polates the pressures of the hadron and quark phases and is 
no more constant. We shall see that the difference between 
the MC and our full calculation causes no serious difference 
regarding the bulk properties of compact stars like the mass- 
radius relation. However, the internal properties are very dif- 
ferent between the two cases, which rnay affect the micro- 
scopic elementary processes in the MP ll36ll . 

If one uses a smaller surface tension parameter cr, the en- 
ergy gets lower and the density range of the MP gets wider. 
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FIG. 5: (Color online) Droplet radius (i?) dependence of the energy 
per baryon for fixed baryon density ps = 0.35 fm~'^ and different 
surface tensions. For any value of R, the density profiles and the cell 
size Rw are optimized. Dots on the curves show the local energy 
minima. The black line shows the energy of the MC case. 



The limit of cr = leads to a bulk application of the Gibbs 
conditions without the Coulomb and surface effects, i.e., the 
so-called Glendenning construction |3]. On the other hand, 
using a larger value of cr, the geometrical structures increase 
in size and the EOS gets closer to that of the MC case. Above 
a critical value of dcrit ~ 70 MeV/fm'^ the structure of the 
MP becomes mechanically unstable i 1 8il : for a fixed volume 
fraction {R/RwY 'he optimal values of R and Rw go to in- 
finity and local charge neutrality is recovered in the MP, where 
the energy density equals that of the MC (see Fig.|5]). 

This mechanical instability is a charge screening effect: 
The optimal values of R and Rw are basically determined 
by the balance between the Coulomb energy (~ R^ in the 
3D case) and the surface energy (~ R^^). However, if 
the charge screening is taken into account, the contribution 
of the screened Coulomb potential is strongly reduced 
when i?, Rw ^ oo. A careful analysis by Voskresensky et 
al. showed that the Coulomb energy changes its sign and be- 
haves like R^^ as i? ^ oo due to the charge screening ifTsll . 
Thus the surface and the Coulomb energy give a local en- 
ergy minimum below (Tcrit ~ 70 MeV/fm^, which disappears 
when the surface energy becomes greater than the Coulomb 
energy above (icrit ~ 70 MeV/fm^. This is in contrast to the 
work of Heiselberg et al. flTt], neglecting the charge screening 
effect, where there is always a local energy minimum at finite 
R. The importance of the charge screening effect has also 
been discussed in the context of the stability of strangelets 

il. 

If we assume uniform density profiles with a given vol- 
ume fraction fv = [R/RwY '^he difference Spch = 

^iPch^ ^ Pch) °f '■^^ charge density between the two phases, it 
is easy to see how the optimal size of the lump is determined 
from the competition between the Coulomb and surface en- 
ergy contributions jsi lssll : The Coulomb interaction energy is 
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in this case 



Vw 



d + 2 



-lid' 



d-2 



(27) 



while the surface energy is simply Es/Vw — dafv/R- 
Hence there is always one energy minimum at finite R for 
a given fy as a consequence of the balance between the 
Coulomb energy and the surface energy 1 17]. 

However, when the Coulomb screening is present, the 
charge density is no more uniform, see Fig. |2] but is rear- 
ranged to attain smaller Coulomb interaction energy. Conse- 
quently the Coulomb potential becomes short range due to the 
Debye screening of the charged particles. The contribution of 
this many-body effect to the energy is then twofold: one is 
the direct contribution of the Coulomb interaction energy and 
the other indirectly via the rearrangement effect of the charge 
density, which is called the correlation energy in Ref. 1 1 8il . 
Considering Pair) as an implicit function of 0(r), the kinetic 
and strong interaction energies can be expanded with respect 
to the particle densities around their uniform values. 



E 



+ 



Vw 



Where £;,Vr = ^^?f + ^^^ 



(28) 

and p° are the quantities for the 
system with uniform densities in the absence of the Coulomb 
screening. Thus the correlation energy gives rise to a new 
i?-dependence in the energy density. The analysis of Voskre- 
sensky et al. showed that the contributions from the screened 
Coulomb potential and the correlation energy exhibit a R~^ 
dependence for large R, and that both have different signs 
fllSn. When the surface energy is added, we can easily see that 
the energy density has a local minimum at finite R as long 
as a is not too large. Actually this local minimum disappears 
above a critical value Ccrit ~ 70 MeV/fm^, which implies a 
mechanical instability of the structure in the MP. 



C. Hyperons in the mixed phase 

One notes in Fig. |2]that no hyperons appear in the MP al- 
though the mean baryon density ps = 0.4 fm^'^ is higher 
than the threshold density for hyperons in pure nucleon mat- 
ter, pb ~ 0.34 fm^'^ (see the black dot in Fig.[T]l. In fact in the 
MC a small fraction of hyperons is present, as displayed 
in the lower panel of Fig. |2l 

This is confirmed in Fig. |6] where we compare the parti- 
cle fractions as a function of baryon density in the full cal- 
culation (upper panel), the MC (middle panel), and the Bulk 
Gibbs construction (lower panel). One can see that the com- 
positions are very different in the three cases, the MP by the 
full calculation lying in between the extreme cases of Bulk 
Gibbs and MC. In particular, a relevant hyperon (E~) frac- 
tion is only present in the hadronic component of the MC. 
In this case the phase transition occurs at constant pressure 
between the pure hadron and quark phases with the equi- 
librium densities pn — 0-34 fm^'^ and pQ = 0.50 fm^'^. 
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FIG. 6: (Color online) Particle fractions of quark and hadron species 
in the MP by the full calculation (upper panel), the MC (middle 
panel), and the Bulk Gibbs construction (lower panel). In the MC 
the phase transition occurs between the pure phases with pH =0.34 



fm~^ and pg = 0.50 fm-^ 



and for a given p b the volume fraction of the quark phase is 
< fv = {pB - Ph)/{pq ~ Ph) < 1- Accordingly the 
hyperon number fraction is always finite in the MC, but in the 
full calculation and the Bulk Gibbs case hyperons are com- 
pletely suppressed. 

This hyperon suppression is due to the absence of the 
charge-neutrality condition in each phase. In a charge-neutral 
hadron phase, hyperons (S^) appear at pB — 0.34 fm^'^. 
This is to reduce the electron Fermi energy by replacing the 
negative charge of electrons with that of particles. In 
the MP, on the other hand, the hadron phase can be posi- 
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FIG. 8: (Color online) Neutron star mass-radius relations for differ- 
ent EOS and three different hadron-quark phase transition construc- 
tions. For the hybrid stars (blue and black curves), the dashed lines 
indicate the Maxwell (upper curves) or bulk Gibbs (lower curves) 
constructions and the solid lines the mixed phase by the full calcula- 
tion. 



FIG. 7: (Color online) Upper panel: Particle fractions of neutral 
matter with electrons (corresponding to neutron star matter). Lower 
panel: The same quantity for charged matter without electrons, the 
low-density part of which corresponds to symmetric nuclear matter. 
Both cases require chemical-equilibrium condition. 



tively charged due to the presence of the negatively charged 
quark phase. This causes the reduction of electron density and 
chemical potential in the hadron phase and, consequently, 
is suppressed. In other words, with the charge-neutrality con- 
dition hyperons appear at low density to reduce the energies 
of electrons and neutrons in spite of large hyperon masses (see 
Fig.|7]l. Without charge-neutrality condition, hyperons appear 
at higher density due to their large masses. 

Thus we conclude that due to the relatively small magni- 
tudes of the surface and Coulomb energies, the EOS of the MP 
is similar to the MC one, but the internal structure of the MP 
is very different. In particular the role of hyperons is strongly 
reduced when we consider the deconfinement transition in hy- 
peronic matter. Above a critical value of the surface tension 
parameter, however, the MC is effectively recovered as the 
physical phase transition. These results should be important 
for physical processes like neutrino propagation and baryonic 
superfluidity, besides the maximum mass problem, which will 
be discussed now. 



V. HYBRID STAR STRUCTURE 

Knowing the EOS comprising hadronic, mixed, and quark 
phase in the form P{e), the equilibrium configurations of 
static NS are obtained in the standard way by solving the 
Tolman-Oppenheimer-Volkoff (TOV) equations yj for the 
pressure P{r) and the enclosed mass m(r). 



dP 
dr 
dm 
dr 



Gme (1 + P/e) (l -h Anr^P/r 



1 - 2Gm/r 



47rr^e 



(29) 
(30) 



being G the gravitational constant. Starting with a central 
mass density e(r = 0) = ec, one integrates out until the 
surface density equals the one of iron. This gives the stel- 
lar radius R and its gravitational mass M = m{R). For the 
description of the NS crust, we have joined the hadronic EOS 
with the ones by Negele and Vautherin Ii39il in the medium- 
density regime, and the ones by Feynman-Metropolis-Teller 
[40] and Baym-Pethick-Sutherland \4V\ for the outer crust. 

Fig.[8]compares the mass-radius relations obtained with the 
different models. The purely nucleonic EOS (green curve) 
yields a maximum NS mass of about 1.82 M©, which is re- 
duced to 1.32 A/0 when allowing for the presence of hyperons 
(red curve). This feature has been shown to be fairly indepen- 
dent of the nucleonic and hyperonic EOS that are used 
The canonical NS with mass of about 1.4 Mq can therefore 
not be purely hadronic stars in our approach. In fact, the inclu- 
sion of quark matter augments the maximum mass of hybrid 
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FIG. 9: (Color online) Internal structure of a 1.4 M© neutron star obtained with three different phase transition constructions. The upper 
panels show total energy density and pressure and the lower panels the overall particle fractions as functions of the radial coordinate of the 
star, using the bulk Gibbs (left panel), the mixed phase by the full calculation (central panel), and the Maxwell construction (right panel). 



Stars to about 1.5 M©. 

More precisely, we compare in the figure results obtained 
with the two quark EOS S 100 MeV/fm'"^ (blue curves) 
and Bcs (pb) (black curves), and involving the different phase 
transition constructions Bulk, Mixed, and Maxwell. In gen- 
eral, the Maxwell construction leads to a kink in the M{R) 
relation, because the transition from a hadronic to a hybrid 
star occurs suddenly, involving a discontinuous increase of 
the central density when the quark phase starts in the core of 
the star The Bulk construction yields smooth mass-radius re- 
lations involving a continuous transition from a hadronic to a 
hybrid star beginning at rather low central density correspond- 
ing to very low NS mass. 

The MP construction by the full calculation lies between the 
two extreme cases, and with our choice of cr = 40 McV/ fm^ it 
is rather close to the Maxwell construction, smoothing out the 
kink of the hadron-hybrid star transition. This transition oc- 
curs generally at a fairly low NS mass, even below the natural 
minimum mass limit due to the formation via a protoneutron 
star |42] and is thus an unobservable feature. 

On the contrary, the maximum mass is hardly affected by 
the type of phase transition: For the Bcs{pb) model the max- 
imum mass is 1.52 Mq, practically independent of the kind of 
phase transition, whereas for the B — 100 MeV/fm"^ model 
there is a sHght variation of M = 1.45, 1.45, 1.50 Mq for the 
Maxwell, mixed, and bulk construction, respectively. 

Whereas the maximum masses are practically independent 
of the phase transition construction, there are evidently large 
differences for the internal composition of the star. This is il- 
lustrated in Fig.|9l which shows the total energy density, pres- 
sure, and particle fractions as a function of the radial coor- 
dinate for a 1.4 A/q NS. One observes with the bulk Gibbs 
construction (left panel) a coexistence of hadrons and quarks 
throughout the whole interior of the star, whereas with the 
MC (right panel) an abrupt transition involving a discontinu- 
ous jump of energy and baryon density occurs at a distance 



r ~ 7.4 km from the center of the star The small contami- 
nation with Yj^ hyperons in the hadronic phase is not visible 
on the scale chosen. The MP with the full calculation (central 
panel) lies between the two extreme cases, hadrons and quarks 
coexisting in the range r « 5.7 - 8.0 km. In both latter cases 
the pure quark matter core has a higher central pressure and 
baryon number and energy densities than the mixed core of 
the first case. 



VI. SUMMARY AND CONCLUDING REMARKS 

In this article we have studied the properties of the mixed 
phase in the quark deconfinement transition in hyperonic mat- 
ter, and their influence on compact star structure. The hyper- 
onic EOS given by the BHF approach with realistic hadronic 
interactions is so soft that the transition density becomes very 
low if one uses the MIT bag model for the quark EOS. 

The hyperon-quark mixed phase was consistently treated 
with the basic thermodynamical requirement due to the Gibbs 
conditions. We have seen that the resultant EOS is little dif- 
ferent from the one given by the Maxwell construction. This 
is because the finite-size effects, the surface tension, and the 
Coulomb interaction tend to diminish the available density re- 
gion through the mechanical instability, as has also been sug- 
gested in previous articles ifTslfTotl . 

For the bulk properties of compact stars, such as mass or 
radius, our EOS gives similar results as those given by the 
Maxwell construction. The maximum mass of a hybrid star is 
around 1.5 Mq, larger than that of the purely hyperonic star, 
w 1.3 Mq. Hence we may conclude that a hybrid star is still 
consistent with the canonical NS mass of IAMq, while the 
masses of purely hyperonic stars lie below it. 

On the other hand, the internal structure of the mixed phase 
is very different; e.g., the charge density as well as the baryon 
number density are nonuniform in the mixed phase. We have 
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also seen that the hyperon number fraction is suppressed in 
the mixed phase due to the relaxation of the charge-neutrality 
condition, while it is always finite in the Maxwell construc- 
tion. This has important consequences for the elementary pro- 
cesses inside compact stars. For example, coherent scattering 
of neutrinos off lumps in the mixed phase may enhance the 
neutrino opacity (3^. Also, the absence of hyperons prevents 
a fast cooling mechanism by way of the hyperon Urea pro- 
cesses |43]. These results directly modify the thermal evolu- 
tion of compact stars. 

Although we have considered the phase transition at zero 
temperature, our study can be easily extended to finite temper- 
ature, which is relevant to protoneutron stars and supernovae. 

In this article we have not included hyperon-hyperon inter- 
actions and three-body forces among hyperons and nucleons, 
since there are still many theoretical and experimental ambi- 
guities. However, some works hav^suggested their relevance 
for the maximum mass problem |44]: if the hadronic EOS is 
sufficiently stiffened by repulsive interactions, the maximum 
mass problem may be resolved. Even in this case, however, 
the quark deconfinement transition may occur and the proper- 
ties of the mixed phase deserve further investigation. 

Finally, we have considered here a very simple quark mat- 
ter model based on the MIT bag model, but there are many 
works about the properties of high-density QCD. Since color 
superconductivity Il34tl or magnetism 14511 in quark matter are 



closely related to the thermal and magnetic evolutions of com- 
pact stars, it should be interesting to take into account these 
effects in the quark EOS for a more realistic description of 
the mixed phase. For example, one may expect 2SC in the 
quark phase, as inferred from Fig. |2l the number densities of 
u and s quarks become similar in the mixed phase, while the 
quark densities are well different in the uniform quark matter 
134, 46, 47]. It would be an interesting possibility and worth 
studying in detail, but lies outside the scope of the present pa- 
per 
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